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Biharmonic  functions  are  much  used  in  the  theory  of  elastic  solids 


A problem  of  long  standing  has  been  to  produce  a formula  involving  only 


the  boundary  conditions  which  gives  a bound  for  the  value  of  a biharmonic 


function  in  the  interior  of  a region.  Using  results  of  Miranda,  this 


AMS  (MOS)  Subject  Classifications:  31A30,  35B30,  73.41 

Key  Words:  Biharmonic  functions.  Dependence  of  solutions  on  boundary 
data,  Approximations 

Work  Unit  Number  3 (Applications  of  Mathematics) 


Sponsored  by  the  United  States  Army  under  Contract  No.lDAAG29-75-C 


MAJORIZATION  FORMULAS  FOR  A BIHARMONIC  FUNCTION 


OF  TWO  VARIABLES 
J.  Barkley  Rosser 


1 . BACKGROUND . We  propose  to  extend  majorization  formulas  given  by 
Miranda  in  [2].  He  was  considering  a biharmonic  function  u(x,y);  that  is. 


(1.1) 


V2V2u 


in  a region  T.  He  allowed  T to  be  multiply  connected,  but  required  the 
curvature  of  the  boundary  (or  boundaries)  to  be  continuous  around  the  boundary, 
FT,  and  the  tangent  to  be  continuously  turning.  On  the  boundary,  values  of  u 
were  specified, 

(1.2)  u = f , 

and  values  of  the  inward  normal  derivative  were  specified, 

du 


(1.3) 


dv 


= g 


He  considered  f and  g as  functions  of  the  arc  length  around  the  boundary. 

He  also  assumed  that  u and  its  first  partial  derivatives  are  continuous  in  T, 
up  to  and  including  (one-sided)  continuity  at  the  boundary. 

His  Theorem  n says  that  if  f = 0,  then  in  T 

d*4)  |u(x,y)  | <_  /2$(x,y)  max|g|  , 

where  the  maximum  of  |g|  is  taken  over  the  boundary  and  $ (x,y)  is  the  func- 
tion which  satisfies 


(1.5) 

inside  T and  has  ^ 


-1 


0 on  FT. 
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Upon  relaxing  the  requirement  that  u be  zero  on  the  boundary,  Miranda 
also  had  a majorization  formula.  He  assumed  continuity  of  f'  and  g,  and 
concluded  that  there  are  nonnegative  constants  K and  K_,  depending  solely 
on  the  region  T,  such. that  in  T one  has  1 2 

(1-6)  |u (x,y)  | <_  Kj  6fmax|g|  +max|f'|]  + (1  + K26)max|f|  , 


where 


is  the  (minimum)  distance  from  (x,y)  to  FT.  See  his  Theorem  VI. 


sizes  of 


This  has  the  shortcoming  that  no  clue  is  available  as  to  what  might  be  the 


Professor  L.  Collatz  has  pointed  out  in  conversation  that,  if  one  is  given 
boundary  conditions  for  u,  it  may  be  possible  to  contrive  a specifically 
given  u whose  boundary  conditions  are  not  greatly  different  from  those  of  u. 
Then,  if  one  had  specific  constants  in  (1.6),  one  could  bound  the  difference 
between  u and  u by  means  of  (1.6).  We  undertake  to  find  in  two  cases 
majorization  formulas  involving  specific  constants  that  can  be  used  as 
Professor  collatz  suggests. 

We  will  consider  the  two  cases  where  T is  a circle  and  where  T is  a 
rectangle,  and  will  supply  majorization  formulas  with  specific  constants;  for 
the  rectangle  we  have  to  assume  in  addition  that  f"  is  of  bounded  variation, 
and  will  find  max|f"|  appearing  in  the  corresponding  majorization  formula. 

We  also  need  a slight  additional  hypothesis  on  u at  the  four  corners  of  the 
rectangle. 


In  Section  2,  we  collect  some  auxiliary  formulas.  In  Sections  3 and  4 
respectively,  we  treat  the  circle  and  rectangle.  The  relevant  majorization 
formulas  are  given  near  the  ends  of  the  two  sections  in  Theorems  1 and  2.  In 
Section  5 we  discuss  the  possibility  of  weakening  the  hypotheses  of  Theorem  2. 

2 . AUXILIARY  FORMULAS . We  collect  here  various  pieces  of  information 
that  will  be  useful  in  subsequent  sections. 

10  | | 

Let  z * re  be  a complex  variable.  If  z <1,  then 


1 

1 - z 


1 + z + z 


Integrating  gives 


2 3 

fn(l  - z)  = -z  - ^ ~ 


Taking  real  and  imaginary  parts  gives 


~ £n(l  - 2rcos0  + r2)  = 


rncosn0 


arctan 


rsinO 
1 - rcos0 


rnsinn0 


2 


■ SS-, 

f.  Wr 

m: 


Taking  r = 1 in  (2.3)  gives 


j.  An  2 + i An(l  - cos0)  = - £ 


cosn6 


/ An(l  - cos6)d6  » - ^ An  2 - 2C  , 
0 2 


where  C is  Catalan's  constant. 


(2.7)  C = l (-1) k (2k;  + l)"2  . 

k=0 

Properties  of  C are  discussed  on  p.  807  of  Abramowitz  and  Stegun  [1],  and  a 
value  to  18  decimals  is  given  on  p.  812.  Rounded  off  to  5 decimals,  it  is 


Prom  (2.5) , we  get 


C a 0.91597. 


/ An(l  - cos0)d0  = 2C  - j An  2 

IT  2 


Taking  r = 1 in  (2.4)  gives 


(2.10) 


tt  o_  _ r sinr 

2 " 2 L.  n 
n=l 


(0  < 0 < 2tt)  . 


Taking  r = -1  in  (2.4)  gives 


(2.11) 


0_  r (-1)  n+1sinnO 
2 = n 


(-IT  < 0 < IT) 


By  contour  integration  we  can  show  that  for  0 < r < R 


(2.12) 


(2.13) 


. 2tt  „2  2 

1 c R - r 
3it  J 9 


dt  = 1 


(R2  ~ r2)  r2lX  R - rcos(t  - 0) 

9wR  i a at  = i , 


■MCflfiH 


where 

(2.14)  l2  =>  R2  - 2Rrcos (t  - 0)  + r2  . 


Indeed  (2.12)  is  worked  out  on  pp.  112-113  of  Whittaker  and  Watson  [3] 
as  an  illustrative  example  of  contour  integration. 

3.  FORMULA  FOR  A CIRCLE.  Let  us  take  Miranda's  region  T to  be  a circle. 
Let  u satisfy  (1.1),  (1.2),  and  (1.3).  Let  the  circle  have  radius  R.  Choose 
polar  coordinates  r and  6,  with  the  origin  at  the  center  of  the  circle. 

Then,  by  a formula  of  Lauricella  (4] , we  have 

2tt  2 2 

(3.1)  u(r ,9)  » ^ / f(t)  dt 

2n  0 i2 


£(hLz^!)  /2n  f . (t)  ginct  - *>  dt 

2lfR  0 i2 


+ 


4nR 


2ir 

/ 


0 


g(t) 
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where  l is  the  distance  from  the  point  (r,9)  to  the  point  (R,t);  that  is, 
2 

t is  given  by  (2.14). 


This  formula  is  usually  cited  with  a minus  before  the  last  term  on  the 
right,  because  Lauricella  was  taking  g(0)  to  be  3u/3r,  whereas  we  are 
taking  g(0)  to  be  the  inward  normal  derivative.  Integration  by  parts  in  the 
second  term  on  the  right  of  (3.1)  reduces  (3.1)  to 


(3.2) 


u(r,0) 


(R2  - r2) 
2irR 


f2V  f(t)  R.-.rcosit_-_0)  dt  + R2 


4ivR 


2 2 ir 

b!  9(t) 


This  formula  must  be  used  with  some  caution.  If , at  0 = 0Q,  f ' (0)  has 

a jump  discontinuity  and  g(0)  is  continuous,  then  the  limit  of  the  inward 
normal  derivative  as  one  approaches  the  circumference  along  the  ray  0 = 0Q  is 


(3.3) 


g<e0) 


<V>  - f,(V) 


This  and  other  idiosyncracies  of  (3.2)  are  discussed  in  Picone  [5];  see 
especially  pp.  216-217  and  pp.  257-261. 

If  we  make  the  same  assumptions  that  Miranda  made  to  derive  (1.4)  and 
(1.6),  namely  that  u and  its  first  partial  derivatives  are  continuous  in  T 
up  to  and  including  the  circumference,  that  will  assure  us  that  f'  and  g 
are  continuous.  Then  (3.2)  defines  the  one  and  only  u satisfying  (1.1), 
(1.2),  and  (1.3)  for  the  circle. 
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As  the  integrands  in  (2.12)  and  (2.13)  are  positive,  one  can  immediately 
conclude  from  (3.2)  by  the  mean  value  theorem  that,  for  some  x(r,0)  and  a(r,6), 


(3.4) 


u = f (t (r,0) ) + 


2R 


g(a  (r , 0) ) . 


As  a consequence,  if  g(o)  is  nonpositive  for  all  a,  we  have 

(3.5)  u <_  max  f . 

Similarly,  if  g(a)  is  nonnegative  for  all  a,  one  has 

(3 .6)  u >_  min  f . 

In  Picone  [5] , at  the  bottom  of  p.  216,  these  surprising  conclusions  are 
attributed  to  Miranda,  who  apparently  used  essentially  the  same  proof. 

From  (3.4),  we  immediately  get  our  majorization  formula. 

Theorem  1.  Let  u be  a biharmonic  function  satisfying  (1.1),  (1.2),  and 
(3*3)  in  a circle  of  radius  R,  such  that  u and  its  first  partial  derivatives 
are  continuous  up  to  and  including  the  circumference.  Then  at  a point  r units 
from  the  center, 


(3.7) 


r2  _ 2 

I u | <max|f|  + --2~  - max|g| 


Clearly  this  is  much  superior  to  Miranda's  Theorem  VI  (see  (1.6)).  Also, 
if  f = 0 it  is  appreciably  better  than  Miranda's  Theorem  II  (see  (1.4)). 
Indeed,  for  the  region  under  consideration,  the  $ of  Miranda’s  Theorem  II  is 


(3.8) 


(R2  - r2)/4  . 


This  $ is  biharmonic.  If  we  put  it  for  u in  (3.7),  we  find  that  it  itself 
is  given  as  an  upper  bound  for  itself.  Clearly  this  is  the  best  possible. 
However,  if  we  use  this  4>  for  u in  Miranda's  Theorem  II,  namely  (1.4),  a 
bound  for  u at  the  center  of  the  circle  is  given  which  is  /2  times  as 
great  as  the  actual  value  of  u at  that  point. 


4 . FORMULA  FOR  A RECTANGLE . If  T is  a simply  connected  region  more 
irregular  than  a circle,  one  may  map  it  conformally  into  a circle.  If  the 
region  has  a smooth  enough  boundary,  the  conformality  may  extend  out  to  the 
boundary,  so  that  normal  derivatives  go  into  normal  derivatives.  If  one  has 
a formula  for  the  conformal  mapping,  it  may  be  tractable  enough  that  one  can 
calculate  factors  of  proportionality  for  the  various  derivatives.  Thus  one  can 
sometimes  convert  (3.7)  into  a majorization  formula  for  the  more  general  region. 

When  one  maps  a rectangle  into  a circle,  conformality  certainly  does  not 
extend  out  to  the  boundary  at  the  four  corners.  In  any  case  the  formula  for 
the  transformation  is  much  too  complicated  to  be  of  much  use.  So  we  give  a 
separate  treatment  for  the  case  that  T is  a rectangle. 


« 
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Let  a and  b be  the  lengths  of  the  sides  of  the  rectangle.  Choose 
coordinates  so  that  one  corner  of  the  rectangle  is  at  the  origin,  and  the 
rectangle  extends  a units  along  the  positive  x-axis,  and  b units  along  the 
positive  y-axis. 

Let  u satisfy  (1.1),  (1.2),  and  (1.3).  We  assume  that  u and  its  first 
partial  derivatives  are  continuous  up  to  and  including  the  perimeter.  Also, 
on  each  side  we  assume  that  f"  is  of  bounded  variation  in  the  closed  interval 
consisting  of  that  side. 

At  the  corners,  strange  things  can  happen.  For  one  thing,  a normal  deriva- 
tive is  not  defined  at  a corner,  nor  do  the  normal  derivatives  along  the  two 
sides  have  to  have  the  same  limits  as  one  approaches  a corner.  However, 
continuity  of  3u/9x  as  one  goes  to  the  corner  means  that  the  limit  of  a normal 
derivative  on  a vertical  side  as  one  approaches  an  upper  corner  must  be  either 
f'  or  -f'  along  the  top  at  the  corner.  Whether  it  is  f*  or  -f'  depends 
on  which  direction  is  taken  as  increasing  arc  length,  and  will  be  different  at 
different  corners.  We  assume  further: 

Boundedness  Hypothesis.  For  each  corner  there  is  a neighborhood  of  that 
2 

corner  within  which  7 u is  bounded. 


As  an  illustration  of  a need  for  a boundedness  hypothesis,  we  cite  the 
following  example.  Consider  the  first  quadrant  of  a circle  of  radius  unity  with 
center  at  the  origin.  If  a function  is  harmonic  inside  this  region  and  is  zero 
around  the  boundary,  it  must  be  identically  zero  by  the  maximum  principle.  But 
note  the  harmonic  function 


(4.1) 


2 . -q  sin20 

v = r sin20 — 


It  is  zero  around  the  boundary,  except  for  an  indeterminacy  at  the  origin, 
where  r = 0.  However,  v is  certainly  not  identically  zero.  It  is  because 
of  the  unboundedness  at  the  origin  that  the  usual  maximum  principle  for  harmonic 
functions  fails.  The  reason  we  invoke  our  Boundedness  Hypothesis  is  to  avert 
a similar  difficulty. 


As  we  said,  we  let  u satisfy  (1.1),  (1.2),  and  (1.3).  Choose  u 
the  harmonic  function  inside  the  rectangle  such  that  on  the  perimeter 


to  be 


(4.2)  u = f . 

As  u is  harmonic,  it  is  biharmonic.  So  u - u is  biharmonic,  and  is  zero  on 
the  perimeter.  Applying  Miranda's  Theorem  II  (see  (1.4)),  we  get 


(4.3) 


|u  - u|  <_  /£$  max|g  - g|  , 


where  g is  the  value  of  the  inward  normal  derivative  for  u,  and  $ is  as 
in  (1.4). 

We  will  later  show  that  u and  its  first  partial  derivatives  are 
continuous  up  to  and  including  the  perimeter.  So  u - u satisfies  those  of 
Miranda's  hypotheses.  However,  Miranda  also  assumed  that  his  region  had  a 
boundary  with  continuous  curvature  and  continuously  turning  tangent.  Lacking 
these,  we  proceed  as  follows. 


i — i ’•  - 


# o 

On  p.  99  of  Miranda  [2]  is  proved  the  Lemma  that  uV  u is  continuous  in 

the  region,  and  approaches  zero  on  the  boundary,  if  f = 0.  One  can  easily 

modify  Miranda ' s proof  to  show  that  if  some  segment  of  the  boundary  has  continuous 
curvature  and  a continuously  turning  tangent,  then  within  a closed  portion  of 
that  segment  uV  u approaches  zero  uniformly  as  one  approaches  the  boundary. 

So  in  the  interior  of  each  side  of  the  rectangle,  one  has  (u  r u)V2(u  - u) 

approaching  zero  as  one  approaches  the  perimeter.  In  the  neighborhood  of  a 

corner,^  7^  is  bounded.  But  u is  harmonic,  so  that  V2u  = 0.  Hence 
9 (u  - u)  is  bounded.  But  u - u approaches  zero,  continuously.  So  we 
conclude  that  (u  - u)V2(u  - u)  approaches  zero  as  one  approaches  a corner. 

So  Miranda's  Lemma  holds  for  the  rectangle,  and  the  rest  of  the  proof 
proceeds  just  as  in  Miranda  [2] . 

For  the  $ of  (4.3)  , we  may  start  with 


(4.4) 


This  is  zero  on  the  two  vertical  sides  of  the  rectangle.  Now  add  to  <|>*  a 
harmonic  function  ^ which  is  zero  on  the  two  vertical  sides , and  equal  to 

1 , 

- — x (a  - x) 

* ** 

on^the  top  and  bottom.  Then  $ = <p  + p . By  the  principle  of  the  maximum, 

<J>  will  be  everywhere  nonpositive.  So  the  <J>  of  (4.3)  will  be  bounded  above 
by  (4.4).  It  will  also  be  nonnegative,  because  it  is  zero  on  the  perimeter. 

By  a similar  argument,  <f>  is  bounded  above  by 


***  1 

<p  = - y(b  - y)  . 


We  have 


|u|  < |u  - u|  + |u|  . 
But,  by  the  principle  of  the  maximum, 

| u | <_  max  | f | . 


So,  by  (4.3)  we  get 

(4.5)  | u | £ max | f | + /2Q  (max | g | + max | g | } . 

So  we  wish  to  find  max|g|. 

Although  f*  will  make  random  jumps  at  the  corners,  f will  be  continuous 
at  each  corner,  since  u is  to  be  continuous  up  to  and  including  the  perimeter. 

In  the  sequel,  we  will  use  superscripts  T,  B,  L,  and  R to  signify  the 
top,  bottom,  left,  and  right  sides  of  the  rectangle.  We  will  express  u as 


(4.6) 


= i + r + 11 


In  this  we  choose  ) 

a polynomial 

(4.7) 

j[  = A + Bx  + Cy  + Dxy  , 

with  A,  B,  C,  and 
of  the  four  corners 

D chosen  so  that  £ has  the  same  value  as  u at  each 
of  the  rectangle.  To  accomplish  this,  we  set 

(4.8) 

A = f(0,0)  , 

(4.9) 

_ f (a,0)  - f (0,0) 

B " a 

(4.10) 

_ f(0,b)  - f(0,0) 

C = b 

(4.11) 

f (a.b)  - f (a, 0)  - f (0,b)  + f(0,0) 

D " ab 

recall  that  f(x,y) 
rectangle. 

is  the  value  assumed  by  u around  the  perimeter  of  the 

We  will  first  determine  the  inward  normal  for  u along  the  top  of  the 
rectangle.  A similar  analysis  will  apply  to  each  of  the  other  three  sides. 


We  have 

- -7-  y = -C  - Dx  . 
3y  L 


f(0,b)  - f(0,0) 
b 


f (a.b)  - f (a,0) 
b 

As  the  right  side  of  (4.12)  is  linear,  its  extreme  values  must  be  (4.13)  and 
(4.14).  However,  (4.13)  is  bounded  below  and  above  by  the  minimum  and  maximum 
of  f',  evaluated  for  x * 0.  (Here  we  are  considering  arc  length  as  increas- 
ing counterclockwise  around  the  rectangle,  so  that  f'  evaluated  for  x = 0 is 

- £ f(o.y)  .) 

Similarly,  (4.14)  is  bounded  below  and  above  by  the  minimum  and  maximum  of  -f', 
evaluated  for  x * a. 


(4.12) 

At  x = 0,  this  is 

(4.13) 

and  at  x = a,  this  is 

(4.14) 
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We  take 


(4.15) 


(4.16) 


(4.17) 


F = max | f | , 

F ' = max | f ' | , 
F"  = max | f " | , 


taken  around  the  perimeter  of  the  rectangle.  Then  we  have  just  shown  that 


(4.18) 


--kl\ ■ 


We  take 


(4.19) 


(4.20) 


(4.21) 


(4.22) 


where 


r = i aT- 

j=l  J sij 


. . DttY 
sinh 


a j TTX 

~ — 7“  sin  J — , 
3irb  a 


. sinh 

IB  = l AB  £— 

j-1  3 sinh  2* 

a 


a l ttx 

— r sin  * — 

irb  a 


. sinli  3’(a ~ xl 

l L - I »L ^ — sin  aa 

j=i  3 sinh  as  b 

D 


. . J MX 

_ ® _ sinh  **7— 

f = l AR  A-  sin  *2.  , 

. . 3 . xna  b 

-1—1  J 1 nVl  X 


b ' 


j=l  sinh 


(4.23)  A j = ~ / si"  ^ |f(x,b)  - f (0,b)  - x J-  ^-  f(0>b)|dx  , 

(4.24)  A®  = \ j sin  jf (x,0)  - f (0,0)  - x f(a,0)  ~ f(0,0)|dx  , 

(4.25)  AJ  = f / sin  ^ |f(0,y)  - f (0,0)  - y f-(-°-'b)  ~ f (0>0) jdy  , 

(4.26)  A*  = | / sin  22  |f (a,y)  - f(a,0)  - y f(-,-b)  ~ f(a,0)|dy  . 


(4.26)  A*  = £ / sin 


« 


mmtm 


We  note  that  I is  zero  on  the  sides  and  bottom  of  the  rectangle,  and 
on  the  top  it  equals 


f(x,b)  - f(0,b)  - x 


f (a,b)  - f(0,b) 


which  is  the  value  cf 


u-I 


along  the  top  of  the  rectangle.  Similar  considerations  apply  to  , £ , and  £ 
with  respect  to  other  sides  of  the  rectangle.  So 

f * f ♦ il  + r 


takes  the  same  values  as 


u - l 


around  the  perimeter  of  the  rectangle.  So 


(4.27) 


i+f+f+r+r 


takes  the  same  values  as  u around  the  perimeter  of  the  rectangle.  But  (4.27) 
is  harmonic.  By  the  uniqueness  of  the  solution  of  a harmonic  equation,  we 
conclude  that  (4.6)  holds. 

We  can  integrate  by  parts  in  (4.23),  (4.24),  (4.25),  and  (4.26)  to  get 


(4.28) 


(4.29) 


(4.30) 


(4.31) 


aT  = - -~r  J f " (x,b)  sin  ^ dx  , 

3 jV  0 


B 2d  r ...  . ..  , jirx  , 

A.  r r J f"  (x,0) sin  dx  , 

t . i 2 ' a 

J ] i 0 


AL  = . f 

j ,2  2 { 

J j ir  0 


/ f " (0,y)  sin  dy  , 


AR  y~2  / (a,y)  sin  dy  , 

^ j ir  0 


where  in  (4.28)  and  (4.29)  the  double  primes  indicate  the  second  partial 
derivatives  with  respect  to  x,  and  in  (4.30)  and  (4.31)  the  double  primes 
indicate  the  second  partial  derivatives  with  respect  to  y.  Be  it  recalled 
that  these  second  derivatives  were  assumed  to  be  of  bounded  variation.  By  the 
result  on  p.  172  of  Whittaker  and  Watson  [3J , the  integrals  on  the  right  sides 
of  (4.28),  (4.29),  (4.30),  and  (4.31)  each  decrease  of  the  order  of  j . So 

each  of  |a^|,  |a®|,  |a^| , and  |aR|  goes  to  zero  of  the  order  of  j 3 as  j 

goes  to  infinity. 
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Thus  the  various  series  on  the  right  of  (4.19),  (4.20),  (4.21),  and 
(4.22)  converge  absolutely  and  uniformly  everywhere  in  the  rectangle,  including 
the  perimeter,  since 


sinh 
a 

sinh 

a 


sin 


< 1 , 


etc.  This  assures  the  continuity  of  £T,  £B,  £L,  and  £R.  If  we  take  3/3x 

or  3/3y  of  J* , £B,  £L,  or  £R,  we  will  multiply  terms  by  a constant  times 
j,  and  replace  some  sines  by  cosines,  or  some  sinh’s  by  cosh's;  the  latter  will 
not  make  an  appreciable  difference  for  large  j.  After  multiplication  by  j, 
the  coefficients  will  still  go  to  zero  of  the  order  of  j-^.  This  will  still 
assure  absolute  and  uniform  convergence  everywhere  in  the  rectangle,  including 
perimeter , so  that  the  partial  derivatives  will  be  continuous.  Needless  to 
say, l and  its  first  partial  derivatives  are  continuous.  So,  by  (4.6),  u 
and  its  first  partial  derivatives  are  continuous  up  to  and  including  the  perimeter. 

Along  the  top  of  the  rectangle,  the  inward  normal  for  £*  is 


cosh 


(4.32) 


m . _ won 

i-  yT  = - y An  flT 

3y  1 .L.  a i . . 

3=1  * sinh 


jnb 


a . jirx 
- — sin  v — 
lirb  a 


We  split  the  right  side  of  (4.32)  into 

^1  + ^2  ' 


where 


(4.33) 


(4.34) 


j = ifi  . 

3=1  J 


we  conclude 


L 

j=l 

a exp|^ 

since 

1**1  . 

I,  » Lim 

1 r*l 

T 

2A. 

-A 


. jrrx 

in  = 


sin 


- 1 


.-2 


where 


By  (4.28),  we  have  for  0 r < 1 

j • jirt  jirx  ^ 
rJsin  d — sin  J — 


l 

j=l 


j 


dt  . 


This  gives 


I,-;/ 


By  (2.3) , we  get 


..j  . _ irj(t  - x)  __i  irj  (t  + x) 

o>  r cos  — * <»  r cos  — “ 

" 4 ~ l 


U-l 


j-1 


j 


Vdt  . 


Ir  “ 2u  ^ «"(t.b){tn(l  - 2rcos  1,(1  + x)  + r2) 

- £n(l  - 2rcos  T-t^  + r2)  jdt  . 

As  f"(t,b)  is  of  bounded  variation,  we  easily  justify  taking  the  limit 
as  r -*•  1.  This  gives 


(4.35)  ^ — J *"(t,b)jin(l  - cos  *(t  + x) ) - £n(l  - 


if  (t  - x)  , , 
cos  1 > dt 


We  have 


J f"  (t,b) £n(l  - cos  *(t  * x)  ) dt 


-hf 


ir+(Trx/a) 


2ff  irx/a 

ir+(irx/a) 


I^J 

2tt  irx/a 


f"( x, b)fn(l  - coss)ds 

v TT  / 


£n(l  - coss)  ds 


1 


Combining  with  (4.36)  gives 


(4.37) 


i-  f I < li£*  . silp. 

3y  L I - \ 2 3b| 


Along  the  top  of  the  rectangle,  the  inward  normal  derivative  for  \ is 


(4.38) 


We  have,  of  course. 


3 rB  r 3"  l 3*x 

- t — ) = ) sin  — . 

3y  L • , a . . 3nb  a 

n=l  sinh  d — 

**  A 


-T-—  < 1 . 

lirb  — 


T-  < h l |ab| 

3y  L — b ' 3 1 
j=l 


Then,  by  (4.29) 


-£lB  I V | f " (x,0)  |<3x  . 

dy  it  b j=l  j 0 


(4.39) 


- -MB  < F..  . 

dy  L - 3b 


Along  the  top  of  the  rectangle,  the  inward  normal  derivative  for  £L  is 


(4.40) 


JL  VL  = - Y 21al 


. _ sinh  ^(a-— ■ 

21  aL  * 


<-l)J  . 


By  (4.30),  this  gives  for  0 < x <_  a 


(4.41) 


j=l  sinh 


n(a---^-  (-l)jsin  ^ 
b £_ 

w lia  j 


Temporarily  denote  the  material  in  the  curly  brackets  by  £ . By  (2.11), 
it  is  just  the  value  at  (x,y)  of  the  harmonic  function  in  the  rectangle 
which  equals  -iry/2b  on  the  left  side  and  is  zero  on  the  other  three  sides. 
So,  by  the  principle  of. the  maximum,  is  nonpositive.  However 

r*  ny  (a  - x) 

L 2ab 

is  a harmonic  function  which  is  nonnegative  on  the  perimeter.  So  by  the 
principle  of  the  minimum, 

V*  iry  (a  - x) 

1 2ab 

is  nonnegative.  So  we  conclude  that 

Using  this  in  (4.41)  gives 


(4.42) 


_3_  pL 
3y  L 


< 

— ab 


/ | f" (0,y) |ydy  < b(a  ~ax)F" 


Our  proof  of  this  required  the  assumption  0 < x a to  insure  convergence 
rapid  enough  to  permit  interchange  of  the  order  of  summation  and  integration 
in  going  from  (4.40)  to  (4.41).  But 

3 rL 
3y  L 

is  continuous  for  0 x a.  So  (4.42)  must  hold  also  for  x = 0. 

A similar  argument  will  give 


(4.43) 


3 rR 
“ 3y  1 


< bxF" 
- 2a 


So  along  the  top  of  the  rectangle,  we  have  the  inward  normal  derivative 
bounded  as  follows 


(4.44) 


3 - 

3y  u 


Along  the  bottom,  one  gets  the  same  bound.  On  each  side,  one  gets  the 
same  bound,  except  with  a and  b interchanged. 

So,  for  a bound  on  our  normal  derivative,  we  have  to  use  whichever  is 
larger  of 


(4.45) 


4Ca  2a  b 
2 3b  2 

n 
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T 


I 


(4.46) 


4Cb  2b_  + a 

2 3a  2 

n 


If  we  subtract  (4.46)  from  (4.45),  the  difference  is  seen  to  be 

^{4.*  a*}  . 

Hence,  we  see  that  if  a ^ b,  then  (4.45)  is  greater  than  or  equal  to 

(4.46) .  So  we  get  the  following  results 

Theorem  2.  Let  u be  a biharmonic  function  satisfying  (1.1),  (1.2),  and 
(1.3)  in  a rectangle  having  sides  of  lengths  a and  b,  where  a ^ b.  Let 
u and  its  first  partial  derivatives  be  continuous  up  to  and  including  the 
perimeter.  On  each  side,  let  f"  be  of  bounded  variation  in  the  closed  interval 
consisting  of  that  side.  At  each  corner,  let  there  be  a neighborhood  of  the 
corner  within  which  v u is  bounded.  Then 

(4.47)  I u I < max|f|  + /2$(max|g|  +max|f'|  + + ^1  max  | f " 1 1 , 


where  the  maxima  are  taken  over  the  perimeter,  and 
which  satisfies 


<J>(x,y)  is  the  function 


(4.48) 


V <|>  = -1 


inside  the  rectangle  and  has  $ = 0 on  the  perimeter,  and  C is  given  by  (2.7). 

NOTE.  If  the  rectangle  is  oriented  with  one  corner  at  the  origin,  one  side 
of  length  a along  the  positive  x-axis,  and  another  side  of  length  b along 
the  positive  y-axis,  then  $ will  be  defined  as  in  the  paragraph  beginning 
just  before  (4.4) . So  $ is  nonnegative  and  is  bounded  by 


jmax{x(a  - x)  , y(b  - y) } 


5.  POSSIBLE  WEAKENING  OF  HYPOTHESES.  In  Thm.  2 we  assume  continuity 
of  both  first  derivatives  up  to  and  including  the  perimeter.  Could  we  relax 
this  assumption  just  at  the  four  corners? 

Tto  get  some  feeling  for  this,  consider  the  unit  square,  ABCD  (see  Fig.  1) , 
situated  in  the  first  quadrant  with  A at  the  origin.  Let  us  require  that  u 
be  zero  on  the  perimeter,  and  ask  for  an  inward  normal  of  -1  along  AB  and 
DC  and  of  +1  along  AD  and  BC. 

Lift  the  figure  out  of  the  plane,  and  flip  it  about  the  diagonal  AC;  AB 
goes  up  and  over  to  the  positive  of  AD,  while  AD  goes  down  and  under  to 
the  position  of  AB.  We  now  have  the  same  boundary  conditions  that  we  started 
with.  So,  if  they  determine  u uniquely,  we  must  have  the  same  values  of  u 
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A B 


Figure  1 


as  before.  However,  since  the  figure  has  been  flipped  upside  down,  the  values 
along  AC  have  been  changed  to  their  negatives.  As  they  come  out  the  same  as 
before,  they  must  be  zero.  Similarly,  we  conclude  that  u is  zero  along  BD. 

As  we  have  an  inward  normal  of  -1  along  AB,  the  values  of  u next  to 
AB  must  be  negative.  Likely  u is  negative  all  inside  the  triangle  AEB, 

sloping  down  from  AB  and  up  to  AE  and  EB.  However,  even  if  this  is  not 

the  case,  consider  what  happens  if  one  starts  vertically  from  AB,  at  a 
distance  x from  A with  x < 1/2.  One  starts  off  with  a slope  of  -1,  which 
certainly  takes  one  to  negative  values  of  u.  But  by  the  time  one  gets  up  to 
AE,  u has  got  back  up  to  zero.  So  one  must  encounter  some  place  of  positive 
slope.  So  the  slope  has  gone  from  -1  to  a positive  value  in  a distance  less 
) than  x.  So,  someplace  along  the  way  32u/3y2  must  be  at  least  1/x. 

2 

To  find  out  what  V u is  doing,  we  have  also  to  get  an  idea  of  the 
behavior  of  32u/3x2.  Let  us  go  along  parallel  to  AB,  and  close  to  it,  from 

AE  to  BE.  We  start  with  u = 0 at  AE  and  finish  with  u = 0 at  BE.  If 

we  are  close  enough  to  AB,  u will  be  mostly  negative  in  between.  This 
indicates  that  32u/3x2  will  tend  to  be  positive. 

Thus  it  appears  that  as  we  approach  A in  the  triangle  AEB,  we  will 
encounter  points  where  V2u  is  greater  than  1/x.  Not  only  are  we  violating 
our  Boundedness  Hypothesis,  but  it  appears  possible  that  uV^  is  not 
approaching  zero  as  we  get  close  to  A.  So  Miranda's  Lemma  is  likely  failing. 
This  voids  our  proof  of  Theorem  2. 

This  suggests  that  if  we  admit  discontinuity  of  first  derivatives  at  a 
corner,  we  may  entail  a violation  of  our  Boundedness  Criterion. 

If  we  retain  continuity  of  both  first  derivatives  at  the  corners,  do  we 
really  need  the  Boundedness  Hypothesis?  The  key  result  is  (4.3).  As  u - u 
is  zero  along  the  perimeter,  the  derivative  along  the  perimeter  must  also  be 
zero.  Given  continuity  at  the  corner,  u - u and  both  its  first  derivatives 
must  approach  zero  continuously  as  one  approaches  a corner.  This  does  not 
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seem  to  leave  much  latitude  for  misbehavior  of  u - u.  We  will  conjecture 
that  this  suffices  to  give  (4.3)  (which  is  adequate),  but  do  not  see  at  this 
time  how  to  prove  it. 

Actually,  though  we  do  not  have  a good  idea  of  the  behavior  of  the  u of 
Fig.  1 (assuming  it  exists,  and  is  unique),  the  best  guess  we  can  make  indicates 
that  it  actually  satisfies  (1.4).  So  perhaps  continuity  of  first  derivatives 
at  the  corners  is  not  really  needed.  However,  we  will  not  venture  to  conjecture 
this.  In  view  of  (3.3)  , a more  likely  conjecture  would  be  that  one  could 
prove  something  like  (1.4),  or  its  parallel  (4.3),  but  with  an  extra  term  on 
the  right  involving  the  differences  between  the  limits  at  a corner  of  a first 
derivative  as  one  approaches  the  corner  along  the  two  sides. 
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